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ABSTRACT 

Method A method for measuring the wettability of rocks by low-field nuclear 

magnetic resonance is described which has application for hydrocarbon reservoir 
engineering or development or in civil engineering . 

The method essentially comprises determining determines the water wet pore 

surface and the oil wet pore surface when the sample is saturated with water and oil, by 
measurements of relaxation times (Tl, T2) of the sample placed in a nuclear magnetic 
resonance device, previously brought to various water or oil saturation states, and 
calculating the wettability index by combination of th«-values the relaxation times 
obtained for said surfaces. 

Application : — hydrocarbon — reservoir — engineering — or — development, — er — m — ewti 
e ngin ee ring, e tc. 
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BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 
Field of the Invention 

fOOOn The present invention relates to a method for measuring the wettability of 

rock samples by measuring the relaxation time by low- field nuclear magnetic resonance. 

The method finds applicati ons no ta bly for a na lysis of rocks tak en fr o m— an 
underground formation containing or likely to contain hydrocarbons. 

Knowledge of various parameters, and notably of the wettability of the rocks, is 
also useful notably for enhanced recovery of a formation, by injection of a fluid under 
pressure, and when the fluid (liquid or gas) best suited for effluents displacement is to 
be determine d b y means of preliminary tests. 

The invention also finds applications in civil engineering for formation hydrology 
in order to evaluate the degree of pollution of formations for example, or in the building 
trade for testing building materials notably in order to determine waterproofing 
treatments for example . 

BACKGROUND OF THE INVENTION 

Description of the Prior Art 
[0002] There are well-known methods for determining the wettability of rocks to 

the water and oil that they may contain, comprising carrying out rock drainage cycles, 

irer- that is displacement of fluids intended to decrease the water saturation, followed by 

imbibition, this term relating to a displacement of the fluids allowing te-an increase in 

the water saturation (Sw) of the rock. The capillary pressure Pc at one point is defined 
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as the difference at equilibrium between the pressure P(oil) of the oil and the pressure 
P(water) of the water. This parameter makes sense only if the two fluids are in the 
continuous phase in the porous medium. For a water wet medium, only the positive 
values make sense. On the other hand, when the medium has a mixed wettability, the 
fluids can remain in the continuous phase for the positive as well as for the negative 
capillary pressures (Pc). 

[0003] For an application of this type, a complete capillary pressure measuring 

cycle generally comprises (Fig.l) : 

a) positive primary drainage of an initially 100 % water-saturated sample (curve 1) ; 

b) positive imbibition (curve 2) ; 

c) negative imbibition (curve 3) ; 

d) negative drainage (curve 4) ; and 

e) positive secondary drainage (curve 5). 

r00041 There are various types of devices allowing the curves of Fig.l to be 

drawn. 

[0005] In a device referred to as-«_a "porous plate" » device, notably described 

in patent US U.S. Patent 4,506,542, the porous rock sample containing two fluids in the 
continuous phase is placed in an elongate cell ended at its two opposite ends by 
capillary barriers permeable to a first fluid. This first fluid is injected under pressure 
through the first membrane and the pressure difference between the injection pressure 
and the pressure of the fluid discharged at the other end is measured. The pressures of 



the two fluids and the capillary pressure Pc are constant all along the sample, and the 
saturation is assumed to be uniform. 

[0006] It is also known to carry out progressive-speed centrifugation by means 

of centrifugation devices such as those described for example in patent s FR French 
Patent 2,772,477 corresponding to fUSf-1 Patent 6,185,985} or FR- French Patent 
2,763,690, or in patents EP European Patent 603,040 corresponding to fUS- U.S. Patent 
5,463,894) or FR- French Patent 2,798,734 filed by the applicant Assignee . 

f00071 The sample initially saturated with a first liquid (brine for example) is 

placed (Fig.4£_8) in an elongate vessel or cup A containing a second fluid of different 
density (oil for example). Cup A is fastened to the end of a rotating arm B and a 
centrifugal force is applied thereto so as to study the displacements of the fluids in the 
sample during at least two distinc t stage s steps . During a first drainage-stage step , the 
assembly is then subjected to a centrifugal force applied along the length of the vessel 
so as to exert on it an expulsion force which tends to cause part of the first fluid to flow 
out. The pressure field created by centrifugation is expressed as a function of the density 
p, of the radius R and of the angular velocity w, by the relation : l / 2 w 2 ** . p (Rmax 2 - 
R 2 ), for each fluid. We impose that the The pressure of the two fluids at the sample 
outlet is the same and is equal to zero at the outlet. At the same time, the second fluid 
flows into the sample. The two fluids move in the sample until a position of equilibrium 
is reached where the force due to the capillary pressure in the pores compensates for the 
centrifugal force exerted. A measuring sonde is placed in the cup on the side of the 
sample. The sonde can be of capacitive type for example and comprise a metal rod 
insulated by a ceramic sheath. The capacity between the rod and the conducting fluid 
(brine), which is proportional to the immersed height, is measured. With this measuring 



means system , the measuring accuracy is 1.5 % of the pore volume. The sonde detects 
the position of the interface between the two liquids in the cup and transmits the 
measuring signals to a measuring signal acquisition and control automaton E including 
hydraulic liquid circulation means and an acquisition device. 

[00081 During the re-imbibition stage, the velocity is decreased so as to study 

the return of the initial fluid into the sample. The local saturations measured with this 
type of device are calculated by an inversion programme from the total amount of water 
expelled from the sample. 

[0009] According to another method, referred to as-^a "dynamic" » method, a 

sample is placed in an elongate cell having water-permeable membranes ended at its 
two ends by water permeable membranes . At a first end, oil under pressure is directly 
injected into the enclosure. Water is also injected, but this injection is carried out 
through the membrane and at a lower pressure. At the opposite end, the oil is directly 
discharged whereas the water flows out through the terminal membrane. Adjustment of 
the oil and water injection rates allows the capillary pressure to be the same at the inlet 
and at the outlet of the enclosure, which leads to a uniform saturation that can be 
deduced from the fluids balance. The capillary pressure is obtained for example by 
measuring the difference between the pressure of the oil and of the water at the 
enclosure outlet. Such a method is notably described by Brown H.W. in-« "Capillary 
Pressure Investigations^ », Petroleum Transaction AIME, vol.192, 1951. Examples of 
implementation are for example described in patents EP- European Patent 729,022 
corresponding to ffJS- U.s. Patent 5,698,7724 or EP-A -EP Patent 974,839 corresponding 
to ft^S- U.S. Patent 6,229,3124 filed by the applicant Assignee . 
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[0010] A method referred to as semi-dynamic method is also known, wherein a 

rock sample imbibed with a first fluid is confined in a closed cell, another fluid under 
pressure is injected at a first end of the enclosure and the opposite end is swept by a 
low-pressure fluid circulated by a pump ing means which carries the drained fluid 
outside. The device compris e s m e ans for m e asuring measures the pressure and the 
saturation of the sample, the amount of fluid discharged and the electric resistivity of 
the sample. This method is implemented for example in patent FR French Patent 
2,708,742 corresponding to ft JS- U.S. Patent 5,679,885} filed by the applicant Assignee . 

1*00111 : Fhese -Once drainage and imbibition curves being- are established, it is 

well-known to calculate the wettability indices WI from the surface areas Ad and Ai 
delimited marked by the positive and negative capillary pressure curves, as shown in 

A 

Fig.l, by the relation WI (USBM) = Log-^. . 

SUMMARY OF THE INVENTION 

[0012] The method for measuring the wettability of a porous rock sample in the 

presence of water and oil according to the invention comprises determining the water 
wet pore surface and the oil wet pore surface when the sample is saturated with water 
and oil, and calculating the wettability index by combination of the values obtained for 
said surfaces. 

r00131 The method finds applications notably for analysis of rocks taken from 

an underground formation containing or likely to contain hydrocarbons. 

[0014] Knowledge of various parameters, and notably of the wettability of the 

rocks, is also useful notably for enhanced recovery of a formation, by injection of a 
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fluid under pressure, and when the fluid (liquid or gas) best suited for effluents 
displacement is to be determined by means of preliminary tests. 

fOOIS") The invention also finds applications in civil engineering for formation 

hydrology in order to evaluate the degree of pollution of formations for example, or in 
the building trade for testing building materials notably in order to determine 
waterproofing treatments for example. 

F00161 Determination of the water wet pore surface and of the oil wet pore 

surface when the sample is saturated with water and oil is obtained for example by 
means of measurements of relaxation times of the sample placed in a nuclear magnetic 
resonance device. 

[00171 The wettability index is for example determined by the relation : 

SM W -SM Q SM w 

7 ^ = e ^ i or the relatlon j nmr = lo gio -rrr- 
SM W + SM o Sm o 

where SM W is the water wet pore surface and SM 0 is the oil wet pore surface when the 

porous medium is saturated with water and oil. 

[0018] According to an implementation mode, the wettability index is 

determined by the following operations : 

a) measuring the characteristic relaxation times of the water-saturated sample ; 

b) measuring the characteristic relaxation times of the oil in the sample in the presence 
of water, in a zone close to saturation (Swir) of the sample ; 

c) measuring the characteristic relaxation times of the water in the sample in the 
presence of oil, in a zone close to residual saturation (Sor) ; 
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d) measuring the relaxation times of the sample in a state where its 100 % oil 
saturation point is reached ; and 

e) combining the measurements of the various relaxation times so as to obtain said the 
index. 

f0019l According to a preferred implementation mode, the relaxation times of 

stages- operations a) to c) are determined after subjecting the sample to centrifugation. 

[0020] According to a preferred implementation mode, the relaxation times of 

stage- operation d) are determined after forced displacement of the fluids in the sample 
placed in a containment cell. 

r00211 An oil whose intrinsic relaxation time (T B ) is as great as possible and as 

close as possible to that of water, dodecane for example, is preferably selected. 

[00221 The characteristic relaxation times are for example those corrresponding 

to either the saturation curves maxima, or to mean values of sa4d -the curves. 

[0023] Experience shows that the measurement of I NMR obtained by means of the 

method is as sensitive but that it requires much less time and that it is applicable to a 
large number of samples. 

BRIEF DESCRIPTION OF THE F1GURES DRA WINGS 

[0024] Other features and advantages of the method and of the device according 

to the invention will be clear from reading the description hereafter of an embodiment 
given by way of non limitative example, with reference to the accompanying drawings 
wherein : 
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[0025] Figure Fig. 1 shows examples of capillary pressure curves 

necessary for determination of the conventional wettability index Iu SB m (United States 
Bureau of Mines)ri 

[0026] Figur e s Figs. 2a and 2b respectively show the distribution of 

relaxation times T 2 for a water wet rock saturated with water and oil, and a 
representation of the phase distribution (matrix in hatched lines, water in light grey and 
oil in darker grey) T ; 

f00271 Figures Figs. 3a and 3b respectively show the distribution of 

relaxation times T 2 for a rock of intermediate wettability saturated with water and oil 
(same central part as in Figure 2), and a representation of the phase distribution (matrix 
in hatched lines, water in light grey and oil in darker grey) ? ; 

[0028] Figure Fig. 4 shows, on capillary pressure curves, the saturation 

values where the NMR relaxation measurements are performed?: 

[0029] Figur e Fig. 5 shows an example of distribution of relaxation times 

T 2 at the saturation values given in Figur e Fig. 4 (from the bottom up : 100 % water, 
Swi, Sor and 100 % dodecane)^ 

[0030] Figure Fig. 6 shows the relation between index I NM r and the 

standard index I US bmt1 

r00311 Figur e Fig. 7 diagrammatically shows a containment cell with a 

permanent heating circuit by fluid circulation^ and 

[0032] Figur e — Fig. 8 diagrammatically shows the layout of a 

centrifiigation device used for implementing the method. 
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DETAILED DESCRIPTION 

DETAILED DESCRIPTION OF THE INVENTION 

[0033] It may b e r e mind e d should be noted that the NMR analysis technique 

essentially consists in applying applies to an object to be tested a first static magnetic 
polarization field Bo intended to intended for aligning the initially randomly oriented 
hydrogen protons nuclei in the direction of the field, then a second impulsive magnetic 
field oscillating at the Larmor frequency, perpendicular to the first one, created by coils 
excited by a control signal to carry out a nuclear magnetic resonance experiment. When 
this impulsive field stops, the return of the nuclei to their initial state or relaxation 
generates electromagnetic signals (echoes) which are detected and analyzed. The 
presence of physical parameters of the object is determined from the amplitude 
characteristics of these signals. 

[0034] Wo propose d e fining a A new wettability index is constructed by 

combining values of the water wet pore surface SM W and of the oil wet pore surface SM Q 
when the porous medium is saturated with water and oil. The index can be calculated 
for example by the relation as follows : 

0°) 

NMR SM W +SM Q 
or by the relation : 

/«.=Jog 10 ^ (lb) 
SM Q 

[0035] Quantities SM W and SM D are obtained for example by measuring the 

dominant relaxation time in a low-field nuclear magnetic resonance experiment. 
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[0036] In fact, in such an experiment, the fundamental relation connecting the 

longitudinal T, or transverse T 2 relaxation time to surface S and to the volume of this 
pore is as follows : 

11 S 

= + Pi i — ( 2 ) 

T T ' V 

X l,2 M \B t 2B V 

where T iB .2B is the relaxation time of the fluid outside the porous medium. Basically, this 
relation comes from the fact that the molecules in the pore undergo diffusion motions 
and interact with the porous wall during the measuring time (the typical relaxation time 
is of the order of 100 ms). The surface interactions are designated by coefficient pi, 2 
referred to as surface relaxivity. Be sid es, we have disregarded a A gradient term is 
discarded that is important when measurements are performed at a high magnetic field 
(> 0.1T). Relation 2 is strictly valid for a pore saturated with a single fluid. In general, 
natural porous media have a pore size and therefore ratio S/V distribution. A relaxation 
times distribution is thus generally observed, but this does not modify the method 
described here. 

r00371 When two fluids are present in a pore within the porous medium, the 

same physical diffusion mechanism is valid, but the distribution of the two fluids in 
relation to the pore surface is of crucial importance. For example, when the medium is 
water wet, the water is at the surface and interacts therewith, whereas the oil is inside 
the pore and does not interact with the surface. When the distribution of the relaxation 
times is measured for such a system, the result of Figure Fig. 2a is observed. The water 
has shorter relaxation times than in the case of a 100 % water saturation, because the 
volume of water V has decreased (see equation 1). The oil is characterized by a 
relaxation time that is the same as when the measurement is performed outside the 
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porous medium (T 1b , 2 b) because there is no interaction with the pore surface. Besides, it 
is well-known that interactions at the water/oil interface do not lead to a relaxation that 
is different from that of the oil outside the porous medium, provided that the oil 
considered is light (refined oil such as dodecane). Diagrammatically, the phase 
distribution is shown in Figur e Fig. 2b. When the same porous medium has a different 
wettability and is not clearly water wet, the distribution of the relaxation times is 
modified ( Figure Fig. 3a). Essentially, the characteristic time of the oil will be shorter as 
a result of an interaction of the oil at the pore surface because the water is no longer the 
closer to the surface. The possible distribution of the two fluids is shown in Figure 
Fig. 3b. Such an effect has been obtained using standard procedures allowing to 
reproduce the surface properties of rocks from a petroleum reservoir : cleaning, then 
aging in the reservoir oil at the temperature of the reservoir at the irreducible water 
saturation. 

[0038] We thus general ize The relation (2) is generalized to a two-phase water- 
oil system. For the water, w e have the relationship is : 

-+A.>ir ( 3 ) 



T T 

A \w,2w 1 \Bw,2Bw 



and for the oil 



1 1 SM o ,„X 



1 lo,2o 1 \Bo,2Bo * ' 



[0039] We can thu s determine Thus quantities SM W and SM 0 can be determined 

by measuring relaxation times T! or T 2 in the porous medium, relaxation times T, B or T 2B 
of the fluids outside the porous medium, liquid volumes V G or V w present in the porous 
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medium. The interaction constants p can be determined separately, but we-it will be 
seen that only the ratio of these constants is necessary. 

[0040] If we introduc e r e lations Relations 3 and 4 are introduced in the base 

relation L w e obtain which results in the relationship : 

Sw(- —)-CpSo(- — ) 

T T T IT 

j w bw 1 o bo sc\ 

1 NMR — -I i ii \ J ) 

SW(- ) + CpSo{ ) 

% 1 bw i o 1 bo 

where Sw==Vw/Vp is the water saturation, 
So=Vo/Vp is the oil saturation, and 
Vp the pore volume. 

Cp = is the ratio of the relaxivities for the water wet (w) and oil wet (o) 

Po 

surfaces. To simplify the notation, indices 1 and 2 have been suppressed. 

f00411 The most suitable saturations for measuring quantities SM W and SM C thus 

remain to be determined. For reasons linked with the calculation of the relaxation time 
distributions, we selee*-the irreducible water saturation (Swi, Figure Fig. 4) is selected 
to determine the oil wet surface in the presence of water (SM W ), and the residual oil 
saturation (Sor, Figure 4) to determine the water wet surface in the presence of oil. We 
also select Also, an oil whose intrinsic relaxation time (T B ) is as great as possible and as 
close as possible to wate r is selected . In fact, if T B is too small, the method will be 
limited to porous media whose ratio S/V is high (small pore sizes), or to media whose 
surface relaxivity is high. Dodecane for example is a refined oil which is suited to the 
proposed measurement because its intrinsic relaxation time (T Bo ) is 1 s, close to the 
relaxation time of water (T Bw , about 2.7 s). Crude petroleum oils should generally not be 
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used because their relaxation times are too small and they also have a relaxation time 
distribution which considerably hinders analysis. The practical details of the calculation 
are given below. 

[0042] From the distribution of the relaxation times at the different saturations 

(see— Seethe example of Figur e Fig. 5), the relaxation time corresponding to the 
dominant peak (Tw at Sor, To at Swi) can be easily determined, and these values are 
used in relation 5. It is well-known that determination of the dominant relaxation time is 
robust and depends little on the distribution calculation process. It can be noted that the 
dominant peak corresponds to that of the oil at Swi and of the water at Sor. The 
presence of the second fluid at these different saturations is important from a physical 
point of view, but it disturbs calculation very little. We can s how - It can be shown that 
whatever the relaxation time of the second fluid present in small amount (water at Swi, 
oil at Sor), the relaxation time of the dominant fluid does not fluctuate much. This 
aspect is important for the calculation robustness. Figur e Fig. 5 shows the distribution of 
the transverse relaxation times T 2 . The longitudinal relaxation time T] can also be used, 
but this measurement is less favourable in general because T t > T 2 but Ti B = T 2B for a 
given fluid. 

[0043] In order to determine the surface relaxivity ratio Cp, we use the dominant 

relaxation times T wl0 o and T ol00 respectively at the two saturations Sw=100 % are used 
and So=100 %. The formula used is as follows : 
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f 00441 Figure Fig. 5 illustrates the two measurements performed. Thus, in total, 

the relaxation times distributions have to be measured at four saturations : 1 00 % brine, 
100 % oil and the extreme saturations Swi and Sor. W e can it is possible to therefore 
use the centrifugation technique to reach, starting from Sw=100%, Swi, then Sor. A 
single rotation stage at maximum speed allows these values to be rapidly reached. For 
measurement at So=100%, we u s e a succession of miscible displacements by injecting 
an alcohol, then the refined oil are used . 

Implementation 

[0045] For implementing the method, w e us e for e xampl e an NMR measuring 

device conventionally comprising (Fig. 7) magnets 1 with an air gap in which coils 2 are 
arranged is used . The coils are connected to an electronic box te- which g enerates the 
excitation signals (signals in the radiofrequency range for example) creating an 
oscillating electric field and to acquire for acquiring the response of sample S to the 
excitation signals. A device of this type is used for example in French patent application 
FR-A-2,823,308. 

Practical implementation example 

[0046] The wettability index is obtained by carrying out for example the 

following succession of stages operation with a reservoir sample whose wettability is 
intermediate : 

1. The sample is 100% saturated with reservoir brine and placed in the NMR 
measuring device (as shown in Figs.7-8 for example) to measure the distributions of 
relaxation times T l5 T 2 and to deduce parameter TwlOO of equation 6 ; 
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2. It is then placed in the centrifugation device (as schematized in Fig.+0_8) and 
centrifuged at maximum speed in a vessel filled with oil until the irreducible water 
saturation Swir is reached ; the distributions of relaxation times T ls T 2 are then measured 
in the NMR device so as to deduce T D of equation 5 ; 

3. After being placed again in the centrifugation device, the sample is centrifuged 
at maximum speed until the residual oil saturation Sor is reached ; it is then transferred 
into the NMR device again to measure relaxation times T l3 T 2 and to deduce T w ; 

4. The sample is thereafter placed in a containment cell such as those described in 
the aforementioned patents and a water and oil miscible solvent, then oil (dodecane for 
example) is injected until the 100% oil saturation is reached. This point being reached, 
relaxation times T„ T 2 are again measured in the NMR device and the value of T 0 ioo 
required for determination of C p is deduced. 

r00471 We thus know Thus all the parameters for calculation of wettability 

coefficient Imm b are known . 

Comparison with index I USB m 

[0048] Comparison between wettability index Iusbm and the new index I NM r 

shows a good correlation between these two quantities (Figure 6) for reservoir rocks. 
The measurement of I USB m can therefore be advantageously replaced by a measurement 
of I N mr, which is as sensitive but requires much less time and is applicable to a large 
number of samples. 

[0049] We have described an An implementation is described example where a 

NMR type relaxometry method is used to measure the surface and the volume of the 




pores. This method is however not limitative. Any other analysis method can be used, 
notably cryomicroscopy. 
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